Upper Limit on the Branching Ratio for the Decay 7r° — > vv 
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A sample of kinematically identified K + — > 7r + 7r° decays obtained with the E949 detector was 
used to search for the helicity-suppressed decay ty° — > vv resulting in an upper limit of 2.7 x 10~ 7 
at 90% confidence level. The upper limit is also applicable to ir° decays into unknown weakly 
interacting particles. 

PACS numbers: 13.20.Cz, 14.40.Aq, 14.60.St 



We report on a search for the rare decay 7r° — > vv 
from the E949 experiment Q, 0, at Brookhaven Na- 
tional Laboratory (BNL). The decay is forbidden by an- 
gular momentum conservation if neutrinos arc purely 
masslcss left-handed particles. A finite neutrino mass 
as evidenced by recent oscillation measurements permits 
the decay to occur. If neutrinos, with mass m„ less 
than half of the 7r° mass, couple to the Z° with stan- 
dard weak-interaction strength, the theoretical branch- 
ing ratio for the 7r° — » vv decay is given as Br(ir° — > 

vv) = 3 x 10~ 8 (m„ / 'm^o) 2 \J 1 — 4 {m v / m^a) 2 for a sin- 
gle Dirac-neutrino type . The experimental upper limit 
for the tau neutrino mass (m„ < 18.2 MeV/c 2 Q) im- 
plies that Br(TT° — ► vv) < 5 x 10 -10 ; cosmological con- 
straints on the neutrino masses imply more stringent 
limits. The branching ratio for 7r° — > vv in the case of 
massive Majorana neutrinos is a factor of two larger Q 
than for Dirac neutrinos because the final state parti- 
cles are identical. In addition to 7r° — > vv, this search 



is sensitive to any decays of the form 7r° — nothing" . 
The 7T° — ►" nothing" decay can arise from several different 
physics processes beyond the standard model, including 
7T° — > vv decay induced by helicity-flipping (chirality- 
changing) pseudoscalar interactions 0, E|, 71-0 — * v\V2 
decay where V\ and V2 are neutrinos of different lepton 
flavor, and 7r° decays to other weakly interacting neu- 
tral states. Astrophysical limits on 7r° vv have also 
been adduced from constraints on the cooling of neutron 
stars through the pion-pole mechanism 77 — * tt° — *• vv 
[To| , altho ugh nuclear medium effects make this model- 
dependent |ll| . 

The current upper limit 01 wa s set by the BNL 
E787 experiment with Br(n° -* vv) < 8.3 x 1CT 7 at 
90% confidence level (C.L.) to all possible vv states. 
A flavor specific search for the decay 7r° — > v^v^ was 
performed by the LSND beam-dump ex per iment, with 
Br(ir° -> v u V u ) < 1.6 x 10" 6 (90%C.L.) 0. 

E949 was designed to measure the rare kaon decay 
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K + — » 7T + ^77 01 . In that measurement, the decay 
K + —> ir + ir° (K^) is a major potential background and 
data is analyzed only with tt + momenta above 0, El or 
below ^fj the iVr2 kinematic peak at 205 MeV/c. In the 
7T° — > vv search, we tag a 205-MeV/c ir° in the detector 
by the presence of a tt + in the kinematic peak. The 
ir° — > w candidates are identified as K^z events with no 
activity other than the K + and ir + in the detector. 

An intense beam of 22 GeV/c protons from the Alter- 
nating Gradient Synchrotron of BNL struck a platinum 
target over a 2.2 s interval (spill) every 5.4 s viewed by a 
beam line 01 with two stages of electrostatic mass sep- 
aration. The typical K + beam intensity (with K + :ir + 
ratio of up to 4:1) at the entrance to the E949 detector 
was 1.3 x 10 7 per spill with momentum 710 MeV/c. Af- 
ter K + 's were discriminated from ir + 's by Cerenkov and 
energy-loss counters, they came to rest in a scintillating- 
fiber target at the rate of 3.5 x 10 6 per spill. The time 
of the 7r + that emerged from the target was required to 
be at least 2 ns later than the time of the incoming K + . 
This "delayed coincidence" requirement guaranteed that 
the 7r + originated from a K + decay at rest, not from a 
scattered beam particle. The momenta of the charged 
decay products were measured in a 1 T magnetic field by 
a drift chamber 0] surrounding the target. The kinetic 
energy and range were measured by a cylindrical array of 
plastic scintillators, the range stack (RS), outside of the 
drift chamber. The resolutions (rms) of the 7r + momen- 
tum (P v + ) , energy (E n + ) and range (R„+ ) from K^i were 
1.1%, 2.9% and 2.9%, resp ectively. Waveform digitizers 
operating at 500 MHz |ljj for the RS readout recorded 
the 7r + — * /i + — > e + decay sequence to distinguish pi- 
ons from muons. Photon detectors covered 47r sr solid 
angle to detect any photon or extra particle from K + de- 
cay. A new photon detection device, the barrel veto liner 
(BVL), was introduced just outside the RS and provided 
2.3 radiation lengths to augment the E787 detector |2ol | 
configuration; with the addition of the BVL, a factor of 
three improvement in the 7r° detection inefficiency was 
expected by Monte Carlo (MC) simulations. Additional 
ancillary photon-detection systems [2l| and an improved 
trigger system were also introduced into the E949 de- 
tector. In 2002, the experiment collected Nk = 1.8x 10 12 
kaons at rest in the target in 12 weeks. 

The 7T° — ► vV search started with the identification of 
decays using the n + kinematics ("K^ tag") in the 
events collected by the K + — ► -k + vv trigger [22J. Selec- 
tion criteria (cuts) on the ir + from the monochromatic 
two-body decay were set at 198 < P n + < 212 MeV/c, 
100 < E n+ < 118 MeV and 28 < R w+ < 33 cm, referred 
to as the "signal box". Potential non-J^ backgrounds 
include K + — > [i + v^ (K^) decays and scattered beam 
pions. These were suppressed and their contribution to 
the total background was estimated using techniques sim- 
ilar to the K + — > ir + vv analysis 0. The decays 
were suppressed with measurements of momentum, en- 



TABLE I: The number of background and candidate events. 
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ergy and range as well as with requirements on the obser- 
vation of the 7r + — > /i + — * e + decay sequence. Beam pion 
background was suppressed by the K + /n + separation in 
the Cerenkov and energy-loss counters and by requiring 
the delayed coincidence in the target. Events with two 
beam particles entering the target, which can defeat the 
delayed coincidence requirement, were rejected by look- 
ing for activity in any of the beam counters at the time of 
the kaon decay. The expected numbers of non-K^ back- 
ground events are summarized in Table [I] Ultimately, 
the search for 7r° — > vV was limited by the detection in- 
efficiency for the 7T° decay photons (20-225 MeV) from 
decay. 

The single event sensitivity for the 7r° — > vv branching 
ratio Br is given by 



SES{tt° -» W) = 



1 



N K Br(K w2 )A K ^ 2 ■ C dls C a 
1 1 



N n CdisCa c 



(1) 



where Br^K^) is the branching ratio of the K V 2 decay, 
Ak w2 is the acceptance of the tag, and N^o is the 
number of 7r°'s collected by the tag. A correction 
factor Cdis was introduced to compensate for the loss of 
K-w2 events from the tagged sample due to the misrecon- 
struction of the n + track by overlapping 7's and e ± 's 
from the predominant ir° decays, which do not occur in 



the 



vv events. The factor was obtained from two 



sets of data produced by MC simulations; one was from 
normal decays, and the other was from K V 2 decays 
where the 7r° — > vv decay was forced. The difference in 
the efficiency of the 7r + reconstruction was used to esti- 
mate the correction factor, Cdis = 1-14 ± 0.01. The cor- 
rection factor C aC c takes into account signal losses due 
to accidental activity in coincidence with the 7r° — > vv 
decay. This factor was obtained from the loss observed 
in a pure sample of K^2 decays (after all activity of the 
muons were removed) by imposing the cut for hermetic 
photon detection (HPD). 

The sensitivity to ir° — > vv was maximized by optimiz- 
ing the parameters for the HPD cut in order to achieve 
the greatest rejection against the n° decay products (77, 
e + e~7) while minimizing the acceptance loss (1 — C acc ) 
due to accidentals. The HPD parameters consisted of 
timing windows and energy thresholds of more than 20 
sub-detectors: typically ±10 ns and 1 MeV. A uniformly 
sampled 1/3 portion of the data ("1/3 sample") was used 
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FIG. 1: Effective 7r° rejection (defined as it rejection x Cacc) 
vs acceptance Cacc of the HPD cut as measured on the 2/3 
sample. The saturated curve at 4 x 10 6 indicates the limit of 
the E949 7r° detection efficiency. 



as a training sample exclusively for tuning the parame- 
ters. To avoid bias, this sample was not used for the sig- 
nal search reported below, nor for the background mea- 
surements shown in Tabled After the parameter space 
was explored to set the HPD parameters, the cut was 
imposed on the remaining 2/3 portion of the data ("2/3 
sample", N v o = 3.02 x 10 9 ) for evaluation. The effective 
rejection (defined as '7r° rejection xC acc ') as a function 
of Cacc for various levels of the cut is shown in Figure ^ 
The HPD cut position optimized on the I /3 sample was 
set at a value of C acc = 0.117 ± 0.002{stat) ± 0.003(sys). 
A total of 99 events were observed in the signal box with 
the final HPD cut. 

Some properties of the 99 events observed in the signal 
box are shown in Figures [5] and [3J The decay time dis- 
tribution in Figure |21 is consistent with the known kaon 
lifetime and does not show any evidence of large contam- 
ination by two-beam background, which would tend to 
flatten the distribution of the surviving events. The 7r + 
momentum distribution in Figure [3] does not show evi- 
dence of a significant contamination by beam pions or 
from K^2 decays. 

The primary reasons for failure to detect photons from 
7T° decay in K„2 events are sampling fluctuations in the 
electromagnetic shower of low energy photons around 20 
MeV and photonuclear interactions of high energy pho- 
tons with undetected products, such as neutrons. While 
the effects of electromagnetic interactions can be esti- 
mated well, there are presently very large uncertainties 
associated with the detailed modeling of detection ineffi- 



FIG. 2: Kaon decay-time distribution with various levels of 
the HPD cut described in the text. All the other cuts except 
the offline delayed-coincidence cut were imposed. The dis- 
tribution was not distorted by the HPD cut confirming that 
the sample was dominated by kaon decays. The depletion 
of events near time zero was due to trigger requirements to 
suppress single beam particle backgrounds. Decay-time fits 
were performed for each plot in a time range of [4ns:30ns]; no 
evidence of two-beam background was found. 

ciencies due to photonuclear processes. Therefore, since 
the overall background contribution from it — > 77 de- 
cays in which both photons go undetected is difficult to 
estimate reliably, we treated all 99 observed events as 
7T° — > vV candidates to set an upper limit. Using Poisson 
statistics, the number of signal events was limited to be 
< 113 at 90% C.L. when 99 events were observed. Sub- 
tracting the non-K^ background of approximately three 
events, the 90% C.L. upper limit of the Br(-K° — ► vv) 
was obtained as : 

Br(n°^vv) < (2) 

v ' 3.02 x 10 9 1.14 x 0.117 v ' 

= 2.7 x 10~ 7 (3) 

The result is three times better than the previous best 
result |l2j . The upper limit obtained above is sensitive 
to any hypothetical weakly-interacting particles, whose 
masses are less than half of the 7r° mass; other deca ys o f 
the kind X° -> "nothing" (e.g. 77, K L , 8 0, and B° jH) 
are experimentally more difficult to measure. 
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FIG. 3: The 7T momentum distribution with various 
levels of the HPD cut. All the other cuts except for the signal 
box cuts were imposed. 
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